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Abstract 

In the framework of the first Born approximation, we present a semirelativistic theo- 
retical study of the inelastic excitation (lsi/2 — > 2s 1 / 2 ) of hydrogen atom by electronic 
impact. The incident and scattered electrons are described by a free Dirac spinor and 
the hydrogen atom target is described by the Darwin wave function. Relativistic and 
spin effects are examined in the relativistic regime. A detailed study has been devoted 
to the nonrelativistic regime as well as the moderate relativistic regime. Some aspects 
of this dependence as well as the dynamic behavior of the DCS in the relativistic regime 
have been addressed. 
PACS number(s): 34.80.Dp, 12.20.Ds 



1 Introduction 

The theoretical study of relativistic electron-atom collisions is fundamental to our understand- 
ing of many aspects in plasma physics and astrophysics. The development of electron-atom 
collision studies has also been strongly motivated by the need of data for testing and develop- 
ing suitable theories of the scattering and collision process, and providing a tool for obtaining 
detailed information on the structure of the target atoms and molecules. Many authors have 
studied this process using numerical tools. Thus, Kisielius et al. [1] employed, the R-matrix 
method with nonrelativistic and relativistic approximations for the hydrogen like He + , Fe 25+ 
and U 91+ ions, where the case of transitions Is — > 2s and Is — > 2p as well as those between 
fine structure n = 2 levels was considered. Andersen et al. [2] have applied the semirelativis- 
tic Breit Pauli R-matrix to calculate the electron-impact excitation of the 2 Si/ 2 — > 2 P° 123/2 
resonance transitions in heavy alkali atoms. Payne et al. [3] have studied the electron-impact 
excitation of the 5s — > 5p resonance transition in rubidium by using a semi-relativistic Breit 



"b . manaut@usms . ma 



1 



2 



Pauli R-matrix with pseudo-states (close-coupling) approach. Attaourti et al. jl] have in- 
vestigated the exact analytical relativistic excitation IS1/2 — > IS 1/2 of atomic hydrogen, by 
electron impact in the presence of a laser field. They have found that a simple formal analogy 
links the analytical expressions of the unpolarized differential cross section without laser and 
the unpolarized differential cross section in the presence of a laser field. 

The aim of this contribution is to add some new physical insights and to show that the 
non-relativistic formalism becomes enable to describe particles with hight kinetic energies. 
Before we present the results of our investigation, we first begin by sketching the main steps 
of our treatment. For pedagogical purposes, we begin by the most basic results of our work 
using atomic units (a.u) in which one has (h = m e = e = 1), where m e is the electron mass at 
rest, and which will be used throughout this work. We will also work with the metric tensor 
g 11 " = diag(l, — 1, — 1, — 1) and the Lorentz scalar product which is defined by (a. b) = a^b^. 
The layout of this paper is as follows. We present the necessary formalism of this work in 
section [2,3 and 4], the result and discussion in section 5 and we end by a brief conclusion in 
section 6. 



2 Theory of the inelastic collision lsi/2 — > 2su2 

In this section, we calculate the exact analytical expression of the semirelativistic unpolarized 
DCS for the relativistic excitation of atomic hydrogen by electron impact. The transition 
matrix element for the direct channel (exchange effects are neglected) is given by 



S fi = -i I dt(ip Pf (x 1 )(j)f(x2) I V d I i) Pl {xi)(pi{x2)) 

r _ 

dt I dr-Lif) (t,r 1 )'fi[) Pi (t } T 1 )(<f>f(x2) I V d I 4>i(x 2 )) (1) 



pi 

+00 



where 

1 z 

V d = 2 

ri 2 n 

is the direct interaction potential, ri are the coordinates of the incident and scattered electron, 
Y2 the atomic electron coordinates, ryi — \ i"i — 1*2 | and t\ =\ ri |. The function ip Pi {xi) = 
ip p (t,Yi) = u(p, s) exp(-ip.x) /y2EV is the electron wave function, described by a free Dirac 
spinor normalized to the volume V, and 4>ij(x2) = 0ij(t,r 2 ) are the semirelativistic wave 
functions of the hydrogen atom where the index % and / stand for the initial and final states 
respectively. The semirelativistic wave function of the atomic hydrogen is the Darwin wave 
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function for bound states [5], which is given by : 

<j>i(t, r 2 ) = exp(-i£ 6 (ls 1/2 )t)^i7(r 2 ) (3) 

where £&(lsi/ 2 ) is the binding energy of the ground state of atomic hydrogen and tp^y (r 2 ) is 
given by : 

^\r 2 ) = (1 4 - ^a.V (2) ) M ( ± Vo(r 2 ) (4) 

it represents a quasi relativistic bound state wave function, accurate to first order in Z/c 
in the relativistic corrections (and normalized to the same order), with ipo being the non- 
relativistic bound state hydrogenic function. The spinors u'^ are such that = (1, 0, 0, 0) T 
and u^' = (0, 1,0, 0) T and represent the basic four-component spinors for a particle at rest 
with spin-up and spin-down, respectively. The matrix differential operator a. A is given by : 



a. A 



For the spin up, we have 
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and for the spin down, we have : 
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where 



(7) 



(8) 



N Dl = 2c/VAc 2 + 1 

is a normalization constant lower but very close to 1. Let us mention that the function <pf(t, r 2 ) 
in Eq. (1) is the Darwin wave function for bound states [6], which is also accurate to the order 
Z/c in the relativistic corrections. This is expressed as 0/(t,r 2 ) = exp(— 2£ b (2si/ 2 )t)<^2^' , (r 2 ) 
with £ b {2si/2) as the binding energy of the 2si/ 2 state of atomic hydrogen. 
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for the spin up and 
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for the spin down. The transition matrix element in Eq. (CQ) becomes : 

/+oo n 
dt / dY l dY 2 %l) pj {t,r 1 )^ ^ Pi (t,ri)<p ] } (t,r 2 )<pi(t,r2)V d 
-oo J 

and it is straightforward to get, for the transition amplitude, 



S 



fi 



—I 



-2irH mel (A)5(E f + S(2s 1/2 ) -E % - S(ls 1/2 )) 



(10) 



11' 



(12) 



where A = \pi — pf \ and 7 is the Dirac matrix. Using the standard technique of the QED, 
we find for the unpolarized DCS 



da 



|P/I 1 



dn f I pi I (4vrc 2 ) 2 \ 2 



J2HPf,s f h°u( Pl ,s t )\ 2 \H mel (^)\ 2 (13) 



SiSf 



3 Calculation of the integral part 



The function H ine i(A) is found if one performs the various integrals : 



H inel (A) = / dr ie * Ari /( ri ) 
Jo 



(14) 



3.1 Integral over r2 

The quantity I(ti) is easily evaluated in the following way. We first write the explicit form of 

J(n) = 

I dr 2 </4 s (r 2 ) 

Jo L r i2 n. 



>i s (r 2 j 



(15) 



Next, we develop the quantity r 12 in spherical harmonics as 
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where r > is the greater of r± and r 2 , and r< the lesser of them. The angular coordinates of 
the vectors ri and r 2 are such that : fj. = (#i,<^i) and r 2 = (0 2 , ipz). We use the well known 
integral [7] 



ryi\ _ _ a^t^ 

d« « m e- a " = -^e— £ — — ite(a) > (17) 



At=0 ^ 

then, after some analytic calculations, we get for : 

^) = ^( ? -4) + ^-- 5 d + ^)r 1 (18) 

3.2 Integral over ri 

The integration over ri gives rise to the following formula : 

P + OO A 

H inel (A) = / rfr ie iAri /( ri ) = --=(/! + I 2 + h) (19) 
Jo v2 

the angular integrals are performed by expanding the plane wave e* Ari in spherical harmonics 
as : 

= ^4 7 rz^(Ar 1 )y, m (A)^(r 1 ) (20) 

lm 

with A = pi — pf is the relativistic momentum transfer and A is the angular coordinates of 
the vector A. Then, after some analytic computations, we get for Ji, J 2 and I3 the following 
result : 

4 /" +0 ° , 3 4 1 



4 /" _3 
11 = 27? io dririe " 3 ^ o(Ari) -27cM(3/2)2 + A2) 

fi 1 f+00 9 1 S 

ft - > + ^) jf " *. - g(l + ^) ((3 *V+1V 

It is clear that the situation is different than in elastic collision [4J , since we have no singularity 
in the case (A — > 0) 



4 Calculation of the spinorial part 

The calculation is now reduced to the computation of traces of 7 matrices. This is routinely 
done using Reduce [8]. We consider the unpolarized DCS. Therefore, the various polarization 
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states have the same probability and the actual calculated spinorial part is given by summing 
over the final polarization Sf and averaging aver the initial polariztion s*. Therfore, the 
spinorial part is given by : 

^£KP/,*/hVp^)| 2 = Tr{ 7 (^c + c 2 ) 7 °(^c + c 2 )} 

S iSf 

2EfEi 



= 2c z 



- (Pi-Pf) + c 



(22) 



We must, of course, recover the result in the nonrelativistic limit (7 — > 1), situation of which 
the differential cross section can simply given by : 



(23) 



da 




128 






(|A nr 


2 + l) 6 



with |A r 



Kj — K/| is the nonrelativistic momentum transfer and the momentum vectors 



(Kj, Kj) are related by the following formula : 

K f = (\K i \ 2 -3/4) 1 ' 2 



(24) 
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- Probability density for radial Darwin (1s) 

- Probability density for radial Dirac 1s 1/3 
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Figure 1: Behavior of the probability density for radial Darwin wave function compared with 
that of the Dirac wave function for small distances and for increasing values of the atomic 
charge number. 
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5 Results and discussions 

In presenting our results it is convenient to consider separately those corresponding to non- 
relativistic regime (the relativistic parameter 7 ~ 1) and those related to relativistic one (the 
relativistic parameter 7 ~ 2). Before beginning the discussion of the obtained results, it is 
worthwhile to recall the meaning of some abbreviation that will appear throughout this section. 
The NRDCS stands for the nonrelativistic differential cross section, where nonrelativistic plane 
wave are used to describe the incident and scattered electrons. The SRDCS stands for the 
semirelativistic differential cross section. 

We begin our numerical work, by the study of the dependence of the probability density 
for radial Darwin and Dirac wave functions, on the atomic charge number Z . As long as the 
condition Za <C 1 is verified, the use of Darwin wave function do not have any influence at 
all on the results at least in the first order of perturbation theory. So, the semi-relativistic 
treatment when Z increases may generate large errors but not in the case of this work. In this 
paper, we can not have numerical instabilities since there are none. For the sake of illustration, 




Figure 2: The long-dashed line represents the 
semi-relativistic DCS, the solid line represents 
the corresponding non-relativistic DCS for a 
relativistic parameter (7 = 1.5) as functions of 
the scattering angle 9. 



Figure 3: The solid line represents the semi- 
relativistic DCS, the long-dashed line repre- 
sents the corresponding non-relativistic DCS 
for various values of the relativistic parameter 
(7 = 1.5, 7 = 2 and 7 = 2.5) as functions of 
the scattering angle 6. 
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we give below the behavior of the probability density for radial Darwin wave functions as well 
as that of the exact relativistic Dirac wave functions for different values of Z . As you may 
see, even if it is not noticeable on the figure 1, there are growing discrepancies for Z = 10 
and these become more pronounced when Z = 20. The QED formulation shows that there 
are relativistic and spin effects at the relativistic domain and the non relativistic formulation 
is no longer valid. 

In the relativistic regime, the semirelativistic differential cross section results obtained 
for the Is — > 2s transition in atomic hydrogen by electron impact, are displayed in figures 2 
and 3. In this regime, there are no theoretical models and experimental data for comparison 
as in nonrelativistic regime. In such a situation, it appears from figures [2] and [3] that in the 
limit of high electron kinetic energy, the effects of the additional spin terms and the relativity 
begin to be noticeable and that the non-relativist ic formalism is no longer applicable. Also a 
pick in the vicinity of 9f = 0° is clearly observed. 

The investigation in the nonrelativistic regime were conducted with 7 as a relativistic 
parameter and 9 as a scattering angle. In atomic units, the kinetic energy is related to 




Figure 4: The solid line represents the semi- Figure 5: The variation of the SRDCSs with 

relativistic DCS, the long-dashed line repre- respect to 9, for various kinetic energies, 

sents the corresponding non-relativistic DCS 
for a relativistic parameter 7 = 1.00053 as 
functions of the scattering angle 9. 
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Figure 6: The variation of the differential Is — 2s cross section of e~ — H scattering at 200 eV. 
The dots are the observed values of J. F. Wiliams (1981) ; the solid line represents the semi- 
relativistic approximation and the long-dashed line corresponds to the non-relativistic DCS. 

7 by the following relation : E k = c 2 (7 — 1). Figure 4 shows the dependence of DCS, 
obtained in two models (SRDCS, NRDCS), on scattering angle 9. In this regime, it appears 
clearly that there is no difference between these models. Figure 5 shows the variation of 
the SRDCS with 9 for various energies. It also shows approximatively in the interval [-5, 
5], the SRDCS increases with 7, but decreases elsewhere. Figure 6 presents the observed 
and calculated angular dependence of Is — 2s differential cross section of e~ — H scattering 
at incident energie 200 eV. Results obtained in two approaches semirelativistic and non- 
relativistic approximations are indistinguishable and in good agreement with the experimental 
data provided by J. F. Williams [9]. 



6 Conclusion 



In this paper we have presented the results of a semirelativistic excitation of atomic hydrogen 
by electronic impact. We have used the simple semirelativistic Darwin wave function that al- 
lows to obtain analytical results in an exact and closed form within the framework of the first 
Born approximation. This model gives good results if the condition Z/c <C 1 is fulfilled. We 
have compared our results with previous nonrelativistic results and have found that the agree- 
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ment between the different theoretical approaches is good in the nonrelativistic regime. We 
have also showed that the non-relativistic treatment is no longer reliable for energies higher. 
We hope that we will be able to compare our theoretical results with forthcoming experimental 
data in the relativistic regime. 
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